o 

(N 

< 



Oh 
I 

^ . 



(N 
> 

m 

00 

o 



% 
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We report on the geometric character of the entanglement dynamics of to pairs of qubits evolving 
according to the double Jaynes-Cummings model. We show that the entanglement dynamics for 
the initial states \%Pq) = cosq:|10) + sinajOl) and |0o) = cosq|11) + sinQ:|00) cover 3-dimensional 
surfaces in the diagram dj x dk x Cu , where Cmn stands for the concurrence between qubits m and 
n, varying < a < 7r/2. In the first case projections of the surfaces on a diagram dj x Cki are conies. 
In the second case the curves can be more complex. We also derive inequalities limiting the sum of 
the squares of the concurrence of every bipartition and show that sudden death of entanglement is 
intimately connected to the size of the average radius of the hypersphere. 

PACS numbers: 03.67.Mn, 03.65.Yz, 03.65. Ud 



INTRODUCTION 



The capacity of quantum systems to entangle is per- 
haps the most intriguing aspect of quantum mechanics 
and is a feature that distinguishes classical from quan- 
tum physics. In a seminal work, Einstein, Podolsky, 
and Rosen [l| have brought this property to discussion 
and since then the subject has been investigated. Re- 
cently, pure bipartite interacting quantum systems have 
proven to be a very useful tool to explore entanglement 
dynamics and unveil several of the intriguing proper- 
ties which govern quantum correlations exchange. Ex- 
amples of such properties are, e. g., the sudden (or 
asymptotic) disappearance of entanglement :2,], the so 
called entanglement sudden birth |3j] , control of entan- 
glement dynamics [J] and entanglement distribution Q, 
an important ingredient for quantum computation. Per- 
haps the best known and explored model is the Jaynes- 
Cummings Model (JCM) [q, where several dynamical 
scenarios have been explored both with and without dissi- 
ation. An analogous model, the Tavis-Cummings model 
has also been used for similar purposes where two 
non-interacting atoms initially entangled, both interact- 
ing with the same cavity mode. The result obtained in 
these two contexts have enlightened entanglement disap- 
pearance in finite time .S -IO J. relations between purity, 
energy and entanglement lllll2l|. invariant entanglement 
|13| and general aspects of entanglement dynamics be- 
tween partitions Il4| -|l7 | . In the present work we show 
that the entanglement dynamics of the Double Jaynes- 
Cummings Model (DJCM) Q exhibits remarkable geo- 
metric properties for the two classes of initial states we 
considered. To quantify entanglement, we use the con- 
currence [l8| . The scenario is a pair of initially entan- 
gled non-interacting atoms "A" and "B" and two cavi- 
ties "a" and "6" which interact locally via the JCM. We 
show that, for initial atomic states belonging to the class 



IV'o) = cosajlO) -|- sinajOl), the relations between con- 
currences describe a conic in a diagram Cij x Cki , with 
ij ^ kl ( ij being equal to Aa, Ah, AB, ab, aB and Bh). 
On the other hand, if the initial atomic state belongs 
to the class \4>q) = cosajll) -I- sinajOO), the geometric 
curve is not as simple. However in all cases when a conic 
is found, the eccentricity can be written as a function 
of the absolute value of the average excitations in A, in 



other words: Vo = 



tr{(TfpQ] 



If the initial atomic state 



is iV'o), Vo gives the probability of the excitation being 
found in only one of the two bipartition Aa or Bb. On 
the other hand, if the initial state is \(t>o), Vq does not 
have the same interpretation. It is important to notice 
that Vq is the predictability which according to the com- 
plementarity relation between two qubits proposed in ref. 
|19j is related to the initial concurrence. Moreover, and 
perhaps more interesting is the three dimensional plot of 
CAsy^CAay^CAb- Each "path" ( i.e. for every alpha, the 
time evolution) is over a surface. We find that this geo- 
metric character can be extended for more dimensions. It 
is possible to define a hypersurface over which the concur- 
rence dynamics between every two pairs i and j defines 
a trajectory over or inside this hypersurface. 

The present work is organized as follows: In section HIl 
we present the physical model, the time evolution for the 
two classes of states, Ii/iq) and \(J)q) ; Next, in section Hill 
we determine the entanglement (quantified by concur- 
rence), and construct the diagram Cij TiCki showing that 
whenever a conic is found its eccentricity is related to the 
predictability as defined in [l9|; In the following section 
IIVI we show the existence of an entanglement surface for 
the dynamics of the pairs of concurrences involving the 
same qubit and justify why the curves of the diagrams 
Cij X Cm will be over that surface; In section |V] we find 
an inequality which describes the entanglement dynamics 
of all qubit pairs. 



II. THE PHYSICAL MODEL 
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Consider a composite system of two identical two-level 
atoms ("A" e "S") and two identical cavities ("a" e "6"). 
The atom "A" ( "i?" ) interacts resonantly with the cavity 



"a" ("6"), respectively, via JCM [6] and the evolution of 
the system is governed by the Hamiltonian 



+ ^ . + , hiO A huj r, 

H = hioa^a + fihj¥b -\ af -\ a, 

2 2 



-gia^a^ 



aa^) ^ g{h'^ a^ + ba] 



(1) 



where a^ [h^ ) and a (b) are the creation and annihilation 
operators of the field inside cavity a (6), respectively. The 
matrices a'^_, a\_ and at, are Pauli matrices of the i-th 
atom, with i — A, B. The cavities are resonant with the 
atoms. So, the frequency of the field inside each cavity 
is equal to the frequency of the atomic transition of the 
atoms internal levels. 



atom 



cavity d 



atom 

B 

cavity b 



We can observe that, at a time t immediately after t = 0, 
the state ^ and © will develop entanglement among 
all the partitions. However, for now we will consider the 
entanglement between qubits (A, B, a e b) and their re- 
lations. To do this we will use as entanglement quantifier 
the concurrence llSlI which is defined as 



C = max 



0,VAi 



A, - ^/A: 



(15) 



where Xi are the eigenvalues of the matrix /o(crj, 
'^y)P*{'^y ® '^y) organized in a descending order. 



III. ENTANGLEMENT DYNAMICS IN THE 



DIAGRAM Cr.iiCk 



kl 



We can easily find the state pij of two qubits taking 
a partial trace over the remaining subsystem. We next 
determine all Cy. 



FIG. 1: A schematic figure of the DJCM. In the left (right) par- 
tition there is the atom "A" ( "_B" ) interacting with the cavity "a" 
( "6" ) , respectively, and there is no interaction between the partition 
"j4a" and "Bfe". 

We consider the cavities initially in the vacuum state 
and some entanglement between the atoms. Consider the 
initial state of the system as 



|?Ao) = (cosa|10) +sina|01))«)|00) 



(2) 



Because of the conservation of the number of excitations 
the time evolution can be determined analytically and 
reads 



\il:t) = xi(i)|10)|00)+a;2(i)|01)|00) + 
-Ka;3(i)|00)|10)+X4(i)|00)|01) 



(3) 



The coefficients will be given by the Schroedinger 
equation, J|V't) = H\ipt), plus the boundary conditions 
a;i(0) — cosa, X2{0) — sina, 0:3(0) = and 2:4(0) = 0. 
They are 

xi{t) = cosacos{gt) (4) 

X2{t) = suYacos{gt) (5) 

Xi{t) = — Jcosasin(f/i) (6) 

X4^{t) — — i sin a sin((7f ) (7) 

Consider also the initial state 

\4>q) = (cosa|ll) +sina|00))® |00) . (8) 

Again for any time 

!</)*) = yi(i)|ll)|00)+2/2(t)|00)|00)+y3(i)|10)|01) + 



-2/4(t)|01)|10)+j/5(i)|00)|ll) 



(9) 



where 



yi{t) = 


e-"^*cosacos^{gt) 


(10) 


y2{t) = 


e'"*sina 


(11) 


V3{t) = 


—i e^*"^* cos a sm{gt) cos{gt) 


(12) 


y4(i) - 


—i e"*"^* cos a sm{gt) cos{gt) 


(13) 


y5{t) = 


~e-"^*coaasm'^{gt) . 


(14) 



A. For the initial state It/jq) 

In this case we obtain 

Cab = |sin2a|cos^(5i) (16) 

Cab = |sin2a|sin2(5i) (17) 

CAa = cos^a|sin(2gt)| (18) 

CAb — I sin2asin(5i) cos(5i)| (19) 

CaB — I sin2asin(5i) cos(5i)| (20) 

Csb = sin2a|sin(25i)| . (21) 

We analyze the geometric structure of entanglement dy- 
namics, in a diagram Cij yiCki- To do this, observe that 
we can sum eq.(|16p with eq.(|17|) and we have 



with Cn 



Cab + Cab — Cq (22) 

sin 2a I is the initial concurrence between 







FIG. 2: Graphic of the straight lines CAB^Cab with a = 
7r/4, 7r/6, 7r/8, tt/IO e 7r/12 for the colors black, blue, red, brown 
and green, respectively. 

the atoms AB. We notice that this equation defines a 
straight line in a diagram Cab x Cab (see figure [2]). The 
lines in equation (|22|) . when a G (0, 7r/2), fill the triangle 
formed by the axis Cab , Cab and Cab + Cab = 1 • Now, 
analyzing equations ([T5)) and ([^0]) we notice that they 
satisfy 



C 



Ab 



Cal 



(23) 



this shows a symmetry between the cavity of one of the 
systems and the atom of the other. Also, dividing ([Tl 
by (|2T|l we easily find 



CAa 



cos^ a 
-r^ — (-^B6 
sm a 



Ct 



(24) 



which is a straight line in the diagram C Aa^^CBh- In 



■a 








0.4 0.6 

AB{ab) 

FIG. 5: Graphic of the semi-ellipse Cyia(„i,) xCao with a = 
7r/4, 7r/6, 7r/8, tt/IO e 7r/12 for the colors black, blue, red, brown 
and green, respectively. The slim violet curve is the scmicircum- 
1. 



ference C^^^^,, 



"C'la 



FIG. 3: Graphic of the straight line CAa^CBh with a = 
7r/4, 7r/6, 7r/8, vr/lO e 7r/12 for the colors black, blue, red, brown 
and green, respectively. The slim violet curve is the straight line 
CAa + Csb = 1- 



2. Concurrence between the atoms (or cavities) versus 

concurrence between one of the atoms and the cavity which 

does not contain it: 



the interval < a < tt/2, the lines ([M]) are limited in 
the region between the lines CAa — 0, Csb — and 
CAa + Csb = 1- Observe that the equations ([^ -[M )) 
define a straight line in their respective diagram C^ x C^i ■ 
The line CAa + Csb = 1 is a superior limit. Using the 
same procedure, and some simplifications, we find other 
conies (ellipses, circumferences and straight lines) which 
we organize as follows: 



1. Concurrence between atoms (or cavities) versus 
concurrence between one of the atoms and its cavity: 



a) CAB(ab)^CBb: 



\CAB(ab) — C^/2) 

cm 



^Bb 

■ 4 

sm a 



(25) 







0.4 0.6 

AB(ab) 

FIG. 4: Graphic of the semi-ellipse CAB(ab)^CBb with a = 
7r/4, 7r/6, 7r/8, tt/IO e tt/12 for the colors black, blue, red, brown 
and green, respectively. 



b) CAB{ab)xCAa: 
{CAB{ab) 



Co/2) 



C^ 

^Aa 



Cl/A 



= 1 



(26) 



{CAB(ab) - Cq/2) + {CaB{Ab)) = Cl / A 



(27) 




0.4 0.6 

AB(ab) 

FIG. 6: Graphic of the semicircumference CAB(ab) ^C^BlAb) with 
a = IT /A, 7r/6, tt/8, tt/IO e tt/12 for the colors black, blue, red, 
brown and green, respectively. 



3. Concurrence between one of the atoms and the cavity 

which does not contain it versus concurrence between one of 

the atoms and its cavity: 



a)Ca 



B{Ab) 



xCAa: 



C 



aB{Ab) 



1 sin2a| 

2 cos^ a 



C, 



(28) 



b) CAb(aB) X Csb 



C 



Ab{aB) 



1 sin2Q| 

2 sin^ a 



-C 



Bb 



(29) 




0.4 0.6 

aB(Ab) 

FIG. 7: Graphic of the straight hne CaB(Ab)''^CAa with a = 
7r/4, 7r/6, vr/S, tt/IO c 7r/12 for the colors black, blue, red, brown and 
green, respectively. The slim violet curve is the semi-circumference 







0.4 o.e 

Ab(aB) 

FIG. 8: Graphic of the straight line CAblaB^^Cgi, com a = 
it/4, 7r/6, 7r/8, tt/IO e it/12 for the colors black, blue, red, brown and 
green, respectively. The slim violet curve is the semi-circumference 

(2Csi,-l)2+(2CA6(afl))' = l- 



To interpret the expressions ([251 - 129)) and their re- 
spective figures (|1]-[S|), it becomes instructive to use the 
predictability 



Vo 



tr {ai 



Po 



(30) 



We use the predictability because it is measurable (the 
module of the mean value of an observable), local and 
is intimately related to the concurrence [19j- For po = 
IV'o)(V'o| we have Vq = |cos(2q;)|, and it is clear that 
Cq -I- Pg = 1. Observe that when T^o = the excitation 
will be equally distributed between the partitions Aa and 
Bb, it will not be localized and the initial entanglement is 
maximum between AB. On the other hand, if Po = 1 the 
atoms will not be initially entangled and the information 
if the excitation will be in partition Aa or Bb will not 
be available. However, we can assure that the excitation 
will be in the partition Aa or in the partition Bb. When 
< Po < 1, all we know is that the excitation has a 
larger probability to be in one of the partitions. 

The eccentricity of the ellipses (P5|) and (|26p can be 
written as a function of the predictability 



Wo 

1+^0 



(31) 



We can determine also the distance / of the focus to the 
center of each ellipse. For the ellipse (pS)) the distance of 



the focus /("^ to its center will be 






7^0(1 T^) 



(32) 



„(a) 



where /< is the focus if < a < 7r/4 and /i. is the 



(a) 



focus if 7r/4 < 
focus /(^^ as 



< 7r/2. The ellipse ^ will have the 



^ib)^,lVo{l±Vo) 



(33) 



f(&) 



where /> is the focus if < a < 7r/4 and f]^' is the 
focus if 7r/4 < a < 7r/2, that is, the opposite case of 
((32|) . This happens because the entanglement of the par- 
tition Aa {Bb) is generated by the evolution JCM and not 
by the initial source of entanglement contained in AB. 
The entanglement generated by the JCM depends on the 
"quantity" of excitation that will be shared between the 
respective atom-field. Thus, when < a < 7r/4 the ex- 
citation, in the state represented by ([3]), will be more 
likely to be found in the partition Aa. Then, the entan- 
glement generated by the JCM in the partition Aa will 
be larger than Bb. This is represented in figure [SJ where 
CAa reaches larger values than 0.5 if < a < 7r/4. In 
this case, the entanglement in the partition Bb has values 
below 0.5, as we can observe in figure S) The same anal- 
ysis is valid in the opposite case, where 7r/4 < a < 7r/2. 
On the other hand, if a £ (0, 7r/2), the eccentricity of the 
ellipses (|25)) and ([26]) are identical, as shown in ([3T|) . but 
the focuses /" and f'^ do not have the same value and do 
not necessarily lie in the same axis, except for the case 
a = 7r/4 when we have circumferences in both cases. For 
example, if a = 7r/6 we have /'' — \/3"/°, where f^ {/'') 
is over the horizontal (vertical) axis, respectively. 

The second case, equation (|27p represents semicircum- 
ferences with radius Co/2. When a — tt/A (Cq — 1) the 
limiting curve is obtained. More generally, we can say 
that a curve defined in its respective diagram Cij x Cki is 
always limited by the semicircumference Cf^ 






The sequential cases, represented by equations 
and ([^5)) . are straight lines with angular coefficient de- 
pendent on the initial entanglement. As we did previ- 
ously, we can write the angular coefficient as functions of 
the predictability. Equation (|28p has angular coefficient 
771 (°) given by 



m 



(a) 



l-Vo 
l+Vo 



±1/2 



,(o) 



(34) 



(a) 



where rrv^ is the coefficient when < a < 7r/4 and 771 
is the coefficient if 7r/4 < a < 7r/2. The straight line 
equation ([^^ will have, in turn, angular coefficient 



,W 



l-Vo 
1 + Po 



±1/2 



(35) 



where, as in the previous case. 



,W 



is the coefficient if 



< a < 7r/4 and m> is the coefficient if 7r/4 < a < 7r/2. 
The opposite occurs for (p4)) . This effect is also due to the 
entanglement given by the JCM as we already discussed 
previously in the eUipse equations ([^ and (PB)) . 



B. For the initial state \<j)o) 

Let us now consider the physical system whose initial 
state is given by equation®. After a time interval t, 
the system will be in state ©. In an analogous way 
as before we determine the concurrences of each pair of 
qubits. Then, we have 



Cab 

Cab 
CAa 

CAb 

CaB 
CBb 
1 



max[0, Co cos^ {gt) ~ 7^] 
max [0, Co sin^ (gt) - 7^] 
cos^(q;)| sm{2gt)\ 

max[0, -Col sin(2.gt)| - 7*] 



max 
,2 



:0,-Co|sin(2.gt)|-7t] 



cos (a) I sin (25t) I 



(36) 

(37) 
(38) 

(39) 

(40) 
(41) 



where 7t = ^ cos^(a) sm^{2gt). Observe that if < a < 

7r/4 we have entanglement sudden death [21 or entangle- 
ment sudden birth [3)]. 

In this case, we have some interesting situations due 
to the symmetry of the system. Notice that the partition 
Aa and Bb will have the same value of predictability. 
Thus, the dynamical entanglement supplied by the JCM 
to Aa or Bb is the same. Observe that Cau = Csb- This 
would not be true if the coupling constant of each JCM 
was different. Due to that same symmetry we also have 
CAb = Csa- Those relations define straight lines (like 
the case of equation ([231)) in their respective diagrams. 
The other diagrams C^ xC^/, however, are not so simple. 
The most interesting case, and also the simplest, is the 
diagram Cab x Cab ■ Consider an instant of time when the 
concurrences Cab and Cab are different from zero at the 
same time. We can the write Cab = Co cos^ {gt)—jt and 
Cab — Co sin^ (gt) — 74 . Notice that using simple algebra 
we can write [Co — {Cab + Cab)] / cos^ {a) = sm^{2gt) and 
{Cab ~ CabY /Cq ~ cos^{2gt). Summing both we have 



{Cab — Cab) Co — {Cab + Cab) 



C2 



+ 



cos^(a) 



= 1 



(42) 



This is a parabola with symmetry axis at 45° of the hor- 
izontal axis {Cab)- On this axis the vertex v is local- 
ized at point u< ~ {0, Co — (1 ± 'Pq)I2} and the focus 
/ at /g = {0, Co T "^o}, where, the index < (>) refers 
to < a < 7r/4 (7r/4 < a < 7r/2), respectively. Be- 
cause of the entanglement sudden death in the partitions 
AB and ab whenever < a < 7r/4, there will only be 
a segment of the parabola in the diagram Cab x Cab if 
the vertex v admits positive values on the axis of the 
parabola. On the other hand, when the vertex is the 
origin or admits negative values, we will only have the 
straight line Cab = or Cab = (observe the figure |9] 
for illustration). If a = arctan(l/2) we have v = {0,0} 
and when the entanglement in one of the partitions dis- 
appears the entanglement of another one resurges, as we 
see in figure [TU] Following the same reasoning, it is clear 
that if < a < arctan(l/2) (or arctan(l/2) < a < ■n/2) 
the entanglement in AB disappears before (or after) it 
appears in ab, respectively (this dynamics is depicted in 
figure [TT]) . 



05 




FIG. 9: Graphic of the parabola Cab^CAB with a = 
Stt/IO, 7r/4, 7r/5 e tt/6 for the colors blue, red, brown and green, 
respectively. 
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FIG. 10: In red (blue) the graphic of Cab (C'at): respectively, if 
a = arctan(l/2). 



Using the sum of equations (1361) and (j37|) . squaring 
them and adding to equations ([38|) or (|4T1) squared, we 
get the following ellipse with expression 



{C^ 



ab 



Cab 



C2 



j ^Aa(Bb) ^ ^ 



(43) 



In figure [T^ we show that we will always have a segment 
of the above ellipse (|43)) , because the entanglement in Aa 
does not suddenly disappear. If < a < arctan(l/2) the 
major semi-axis will be parallel to CAa(Bb)- When a = 
arctan(l/2) we have a circumference and if arctan(l/2) < 
a < 7r/2, the major semi-axis will be parallel to Cab — 
Cab- Then, the eccentricity will be 



5Po-3 

(1 + Vo) 



e= < 



3-5-p, 







4(1 - Vo) 
3 -f 57-^0 



, se < a < ao 



, se ao < c^ < 7''/4 



, se 7r/4 < a < 7r/2 



(44) 



L V 4(1 + Vo) ' 
where ao = arctan(l/2). The focus will be 



' v/(57'o-3)(l + 7'o)/2, se < a < ao 



f^{ 



v/(3-5Po)(l + 7^o)/2, se ao < a < ^/4 (45) 



^ ^(3-^57'o)(l-n)/2, se 7r/4 < a < 7r/2 



u 




FIG. 11: In black (blue) the graphic of Cab (Cab) when a = 
7r/9 < arctan(l/2), respectively. In green (red) the graphic Cab 
{Cab) when a = 7r/5 > arctan(l/2), respectively. 




0.4 0.6 

Aa (Bb) 

FIG. 13: Graphic of the parabola CAa(Bb) ^ {Cab + Cab) with 
a = Stt/IO, 7r/6, arctan(l/2) e arctan(l/3) for the colors black, 
red, brown and green, respectively. 
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0.4 0.6 

Aa (Bb) 

FIG. 12: Graphic of the ellipse CAa(Bb) ^ [Cab — Cab) with a = 
7r/4, 7r/6,7r/8, vr/lO e 7r/12 for the colors black, blue, red, brown 
and green, respectively. 



Observe that when a — ao we have Vo = 3/5, / ~ 
and the ehipse becomes a semicircumference. Notice 
that if < a < ao the entanglement of AB disappears 
before the appearance of entanglement in ab. However, 
the entanglement of Aa is given by the JCM and does 
not remain zero in any finite interval of time. As a result, 
there will be a time interval such that Cab ~ Cab will be 
zero but the entanglement between Aa will not. Thus, 
CAa will admit values larger than Cab — Cab and we have 
the major semi-axis parallel to CAa- Consider now the 
expression used previously, [Cq — {Cab + Cab)]/ cos^ (a) = 
sm'^{2gt). Using equation (pS)) or (|4T1) we arrive again to 
a parabola whose equation reads 



[Cab + Cab) — Co — 



^Aa{Bb) 



(46) 



In figure [13] it becomes clear that the vertex v and the 
focus / are localized on the axis [Cab + Cab) at points 
given by w = {0, ^1 - {Vo)^} and /g = {0, VT^^ - 
(1 ± T'o)}- As before, the sub-index < (>) is if < 
a < 7r/4 {tt/A < a < tt/2), respectively. We always 
have a segment of this parabola in the diagram [Cab + 
Cab) xCAa{Bb)i bccause its vertex is limited between 
and 1. We also know that [Cab + Cab) will not be zero 
if ao < a < 7r/2 and in this interval the parabola does 
not touch the axis CAa{Bb) in this cases. 

From equations ([55]) or (|^T|l we can write | sm{2gt)\ = 
CAa{Bb) cos^ a and substitute in ([55]) or (HOI) . With some 



simplifications we have 

CAb(aB) + T; 5 — CAa(Bb) 

^ ' 2 cos^ a V ^ ' 



c^ 

2 






8 cos^ a 



(47) 



that, like the previous case, is also a parabola with vertex 
V and focus / localized at points 



h = 



Vi^^ TT'o 



the sub-index follows 
parabola of equation 



the previous notation. 
47l) will touch twice the 



The 



— CAa{Bb) when < a < ao. This happens because if 
< a < ao, there will be entanglement sudden death in 
the partition CAb(aB)- If ao < a < 7r/2, on the other 
hand, there will not be sudden death and the segment of 
the parabola only touches the axis CAa(Bb) at the origin, 
as showed in the figure [TH 










0.4 0.6 

Aa (Bb) 

FIG. 14: Graphic of the parabola CAa(Bb)^CAb(aB) with a = 
7r/3, Stt/IO, 7r/4, it/6, arctan(l/2) e arctan(l/3) for the colors black, 
red, brown, green and orange respectively. The violet curve is the 
curve II47I I for gt = 7r/4 and varying a of until tt/2. 



IV. THE ENTANGLEMENT SURFACE 



In the previous section we explore the diagram 
CijxCki for two different initial states. Because of the 



unitary evolution of the physical model and the existence 
of an entanglement invariant [13l] , it becomes instructive 
to analyze the three dimensional diagram Cij x Cik x Cu 
for the i-th qubit. First we analyze such diagram for 
the atom A. For the initial state ([2]), the concurrences 
between the atom A and any other qubit are given by 
equations ([T5)) , (IT5|) and (|19l) . If we make the paramet- 
ric graphics of this concurrences we have curves, for a 
determined value of a, in a diagram C^sxCyiaxCAb, 
as showed in the figure [151 Naturally, if we look at 



agram Caj x C^fe will give the curves of that respective 
diagram. As in the previous case we can draw all possible 




FIG. 15: Graphic of the diagram Cab ^CAa^CAb for the atom 
A and initial state l(2j. From the superior curve to the inferior, we 
have, respectively, a = 7r/4, 7r/6, 7r/8, tt/IO c it/12. 

the projections of this curves in the planes C'AB'x^CAa, 
Cab X CAb and CAa x C^b we get the graphics drawn 
in figures [5l [6] and [71 respectively. Now, interestingly 
enough, if we draw all the possible curves (varying a of 
until 7r/2) in the diagram C^s x C^iq x C^b we have a 
surface in that space depicted in figure [THl A point over 





FIG. 17: Graphic of the diagram Cab ^CAa^dAb for the atom 
A and initial state lO, with a = 7r/4, n/6, 7r/8, tt/IO e 7r/12. 

curves in the diagram C'ab x CAa x CAt if we vary a of 
until 7r/2 and we find an entanglement surface, see the 
figure [m Points over this surface also gives us how much 




FIG. 18: Entanglement surface of the atom A in the diagram 
Cab ^dAa xCyK, for the initial state ^. 

entanglement there is in each of the subsystems AB, Aa 
eAb. 

This same conclusions are also true for B, a and b. So, 
in a general way, we can say that any trajectory in the 
diagrams Cij x Cik , Cij x Cu and Cik x Cu will belong to 
the surface in Cij x Cik x Cu and will be projections in its 



respective diagrams, where i, j, 
{A, B, a and b) of the system. 



k and / are the 4 qubits 



FIG. 16: Entanglement surface for the atom A in the diagram 
Cab ^CAa^CAb for th^ initial state |[2ll. 

that surface informs us how much entanglement there is 
in each one of the partitions AB, Aa and Ab. If now 
we consider the initial state ([5]) and draw the paramet- 
ric graphics, for determined values of a, in a diagram 
C/iB xC/ia xC^h, we also have, curves in that diagram, 
as depicted in figure [T71 where its projections in each di- 



HYPERSPHERE SHELL OF THE 

ENTANGLEMENT DYNAMICS 



Next, we are going to use a result already obtained in 
[5| and [ig]. In these references, they observed that for 
the initial state (0 we have C^s + C^{, + C^^ + C^b = C^. 
Without loss of generality, we can sum in both sides the 



term C\^ 



'^Bb *^^^ ^'^i^ §i'^^ u^ ^AB' 



'^Ab' 



'C„T^ + C„h + 



C^^ + C^f^ = Cq + (cos a + sin a) sin {2gt). This last radius of the hypersphere given by 



expression can be transformed in the inequality C\, 






^aB 



C^ab 



^Aa 



^Bb 



< C'o + (cos* a + sin"* a). 
Now, if we use simple trigonometric relations and the 
predictability, we can rewrite this equation as 



^AB + ^Afa + C'o 



(J2 

^ab + ^Aa + ^Bb < 1 + -;r 



(48) 



which is a hypersphere with radius \/{2 + Cq)/2 in a 
space where the axis are the concurrences between pairs 
of qubits. On the other hand, we can generalize the above 
inequality to 



Co < Cj^g+Cj^fj 



-CIb + CI, + C\, + C%, < 1 + ^ (49) 



which defines a limited region (a hypersphere shell) inside 
the hypersphere defined by ea. (l48l) . Thus, any curve in 
a diagram where the axis are concurrences between pairs 
of qubits and the initial state is ([2]) will lie either on 
the surface or in the interior of the hypersphere shell 
((49)) . So, we can speculate that, in the same way that the 
curves in the diagrams Cy- x Cik are projections of curves 
in Cij X Cik X Cii , the surface defined in Cij x Cik x Cn is 
a projection of the surface of a hypersphere that is in a 
space of greater dimension. 

We can make the same analysis for the initial state 
([5]). However, in that case jS, 16] we have only the in- 
equality < Cls + C^b + CIb + Clf, < Ci and, as 
done previously, we can sum both sides with the term 
Cla + C%h = [1 + cos{2a)]^sin^{2gt)/2. With a simple 
algebra we can express the result of this sum in the in- 
equality < C^^ + C\, + Cl^ + Cl, + C\, + C|, < 
Cl + [1 + cos(2a)]V2. We have the predictability Vq 
equals to cos(2q:) if < a < 7r/4 and — cos(2a) if 
7r/4 < a < 7r/2. Using this and Cg 
rewrite the inequality 



T^Q = 1 we can 



< cls 



" ^Ab + Cq 



' C^ab + C^Aa ' 



ck<i 



c^ 



■±Vo 

'. ^^^^ 
where on the right hand side of the equation we will have 

(1 + Cl/2 + Vo) when < a < 7r/4 and (1 + C^/2 - Vo) 

when 7r/4 < a < 7r/2. This inequality must be valid 

during the whole evolution and, in a space defined by 

axis corresponding to the concurrences Cij, we have the 



R^ = 



C^ 



(51) 



where we will have Ry {R<) when < a < 7r/4 (7r/4 < 
a < 7r/2), respectively. It is noteworthy that for < 
a < 7r/4 there will be sudden death of entanglement in 
some partitions. On the other hand, for 7r/4 < a < tt/2 
there will not be sudden death for any partition. Thus, 
we have i?> when there is sudden death and i?< other- 
wise. Note that for < a < arctan(l/2) there will be a 
time interval At = [arccos(-\/tana ) — arcsin(-\/tan a Ml g 
where Cab = Cab ~ CAb = CaB ~ (as observed in Q). 
Since the hypersphere is defined by the concurrences be- 
tween pairs of qubits, one would intuitively expect, in this 
conditions and during the time interval At, to obtain _R< 
in place of i?> since only Cao and Cgb are different from 
zero. The fact that the average radius of the hypersphere 
increases is a consequence of the dynamical entanglement 
CAa and Csb- When < a < 7r/4 the entanglement of 
the partitions Aa and Bh will attain maximum values be- 
tween 1/2 and 1. Thus, the maximum value of C\^ + C^f^ 
willbe between 1/2 and 2, contributing substantially to 
the inequality (JSO]) . 



VI. CONCLUSIONS 

We presented a detailed study of the geometric char- 
acter of the entanglement dynamics of to pairs of qubits 
evolving according to the double Jaynes-Cummings 
model. Although, this is a analytically solvable simple 
model, it exhibits a very rich dynamical structure which 
we explored here in order to give a geometric meaning to 
the entanglement dynamics. As it became clear, its very 
difficult to generalize our results to other more sophisti- 
cated models or initial conditions. However, we strongly 
believe that there is an intimate connection between the 
average radius of the hypersphere and the phenomenon 
of sudden death of entanglement. We hope this work has 
shed some light on the geometry character of quantum 
mechanics. We also hope experimentalists will become 
sufficiently intrigued with the results given that the dou- 
ble Jaynes-Cummings model is within todays available 
technology. 
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